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Monolayers of coUoidally synthesized cobalt-platinum nanoparticles of different diameters char- 
acterized by TEM (transmission electron microscopy) were deposited on structured silicon oxide 
substrates and characterized by SEM (scanning electron microscopy), GISAXS (grazing incidence 
x-ray scattering), and electric transport measurements. The highly ordered nanoparticle films show 
a thermally activated electron hopping between spatially adjacent particles at room temperature 
and Coulomb blockade at low temperatures. We present a novel approach to experimentally deter- 
mine the particles charging energies giving values of 6.7-25.4 meV dependent on the particles size 
and independent of the interparticle distance. These observations are supported by FEM (finite 
element method) calculations showing the self-capacitance to be the determining value which only 
depends on the permittivity constant of the surrounding space and the particles radius. 
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Nanostructured materials have raised substantial in- 
terest in research and industry P, 0, H, 0, They 
are used or discussed for a wide range of apphcations 
0, hke in optical devices Hi, or in electronic Q and 
magnetic architectures [§, For the chemical indus- 
try nanostructured catalysts are of tremendous interest 
[Tol . [Til [H, • Also under intensive research are appli- 
cations in sensors [13, [iH EBl and medicine [13, • Fur- 
thermore, nanoparticles and assembhes of nanoparticles 
are interesting model systems for investigations of the 
matter properties on the nanometer scale [H, [ll| . Size 
and form dependence of physical and chemical proper- 
ties and photophysical processes in constrained systems 
are still under investigation [2^, [111 • While the building 
blocks of conventional macroscopic crystals are atoms or 
molecules, nanoparticles can compose "artificial solids" 
or "supercrystals" [H, [l^l ■ The use of nanoparticles as 
building blocks allows control over certain properties of 
these novel materials. The electrical transport proper- 
ties of these systems are the subject of numerous pub- 
lications [1 [3, [1, [H [H [H, [111. Ordered monolayers 
of metal nanoparticles represent a simple system which 
allows verification of theoretical models. A vast number 
of experimental and theoretical results dealing with the 
transport prop erties of the nanoparticle arrays has been 

published [M IM 111, [M, [Ml, [13, [11, [13, IM IM [13 

(to mention only a few). Despite this tremendous work 
many features of this systems are still unclear. 

In this paper we discuss the different transport prop- 
erties of nanoparticle films and supercrystals which are 
governed in contrast to conventional crystals or metal 
films by the presence of tunnel barriers between the par- 
ticles if the tunneling resistance (i?t) considerably ex- 
ceeds the resistance quantum Rt ^ = 25.8fcil, with 
the Planck constant h and the elementary charge e. At 
low temperatures small nanoparticles show current sup- 
pression known as the Coulomb blockade [11] which orig- 




FIG. 1: High resolution TEM images of the monodisperse 
cobalt-platinum nanoparticles used in this work (upper fig- 
ures). The big particles were cubic in contrast to the other 
nanocrystals. SEM micrograph of a cobalt-platinum nanopar- 
ticle monolayer deposited on an substrate structured with 
gold electrodes (in top and bottom area of the lower picture) . 
The inset shows the obtained Fourier transformed pattern. 



inates from a charging energy Ec of the particles due to 
their very small capacitance. The total capacitance of 
the particles Cs is the sum of the self-capacitance Cs 
and the contact capacitance of a particle to the adjacent 
particles Cgeo- 



2Cs 



with Cs = + Cgco (1) 
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The observation of the Coulomb blockade requires the 
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thermal energy to be lower than the particles charging 
energy: > ^6^. 
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FIG. 2: GISAXS pattern and IsGISAXS simulation of cobalt- 
platinum particles on silicon substrate. Simulation model: 
2D hexagonal lattice (nearest neighbor distance = 8.8 nm) of 
flattened spheroidals on a substrate with lateral diameters 
of 7.6 nm and a height of 6.5 nm. Lower graph: GISAXS 
out of plane scattering curves (at Yoneda maximum) from 
cobalt-platinum particles on silicon. Curves were analyzed 
using a 2D hexagonal lattice function {q = 8.8nm, spherical 
form factor with diameter of 7.6 nm) . The calculation of the 
GISAXS out of plane cuts along qy is based on a model which 
does not take uncorrelated roughness and diffuse scattering 
into account. Directed to low q these properties lead to an 
increasing intensity. 

A bottom-up approach was used to synthesize cobalt- 
platinum nanoparticles by wet chemical methods for- 
merly developed in our group (soj . We managed to syn- 
thesize particles of different diameters by injecting the 
cobalt precursor solution at appropriate temperatures of 



the platinum precursor/ligands solution. The resulting 
nanocrystals showed a chemically disordered fee lattice 
as was shown by XRD and electron diffraction. The com- 
position of Coq. 14-0. 22Pto. 86-0. 78 for different syntheses 
was determined by EDX (energy dispersive x-ray analy- 
sis). The particle diameters were investigated by TEM, 
the obtained values were 4.4 nm, 7.6 nm, and 10.4 nm ([1]). 
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FIG. 3: As an example the I-V curves of the 7.6 nm particles 
at different temperatures are shown. The films consisting of 
other particle sizes showed qualitatively the same behavior. 
The I-V curves showed distinct deviations at different tem- 
peratures. While the curve at room temperature is almost 
linear (upper graph) a clear current blockade is present at 
4.3 K (lower graph). While at room temperature our data 
followed a sinh function, the Nordheim-Fowler equation was 
applicable at low temperatures. The inset at the upper graph 
shows the development of the experimental zero bias conduc- 
tance of the same film with temperature. The corresponding 
fit was calculated using the simple thermally activated model. 
The inset of the lower figure shows the evolution of the differ- 
ential conductance with temperature for the 7.6 nm particles. 

For GISAXS and electrical measurements we deposited 
the nanoparticles from the toluene solution on silicon 
oxide substrates by the Langmuir-Blodgett technique as 
we have demonstrated recently [1^. SEM micrographs 
showed superlattices of hexagonally ordered monolayers 
with a few small cavities (see [T]) . 
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The microscopic methods SEM and AFM (atomic force 
microscopy) are not suitable to investigate the surface 
morphology on the scale of several millimeters. Micro- 
graphs show only a small area of the sample and for 
better statistics additional methods are needed. Here, 
GISAXS is the first choice, as it reveals information from 
a large surface area (several mm^) of the surface. The 
scattering curves can be compared with the Fourier trans- 
formations of the microscopic images, to demonstrate 
that the same nanostructure is uniformly spread over a 
large surface area. From simulations of GISAXS patterns 
it is possible to get information about the form factor 
and the interference function, which leads to the size and 
shape of the particles and their lateral distances [io! ]. 
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FIG. 4: The determined film conductivity scaled exponen- 
tially with the interparticle distance (upper figure). Field 
dependence of the activation energy was well-fit by an expo- 
nentially decaying equation (lower graph). 

The measured GISAXS patterns can be compared 
with simulations done by the IsGISAXS software [4ll |. 
[2] shows such a pattern and a simulation for cobalt- 
platinum nanoparticle monolayer. Assuming the de- 
picted model of flattened spheroidal particles packed in a 
two-dimensional array the simulation leads to an almost 
identical pattern. The dark line indicates were the detec- 
tor was protected with a lead stripe to prevent damage 



due to high scattering intensities. This method allows 
determining the lateral distances between the particles 
and their height and shape, which are identical with the 
calculations made by analyzing the intensity cut along qy 
at the critical angle of the substrate [42]. This analysis 
is shown in[2l GISAXS allowed us to determine not only 
the long range ordering in our films but also to determine 
precisely the shortest surface to surface interparticle dis- 
tances. 

In order to probe the transport properties of the 
nanoparticle films we measured the DC response at differ- 
ent temperatures. Although the big particles were cubic 
no form induced effects were detected in our electrical 
measurements. At room temperature the obtained I-V 
curves are point-symmetric and almost linear. Toward 
lower temperatures, however, the current- voltage charac- 
teristic becomes more and more nonlinear and the device 
is clearly in Coulomb blockade regime at moderate volt- 
ages as can be seen in |3l Due to the nonlinearity of the 
I-V curves at any temperature and the smooth transition 
from the blocked regime into unblocked with rising volt- 
age it is not appropriate to extract the blockade voltage 
using linear interpolation. Also the transition from noise 
level to measurable current levels was not applicable to 
determine the Coulomb blockade as the blockade voltage 
determined this way was dependent on the measurements 
temperature and can be expected to depend on the noise 
level of the used experimental setup. 

At room temperature the I-V curves could be fitted 
well by a sinh function suggested in literature for elec- 
tron hopping between spatially adjacent energy states in 
semiconductors [43|. In contrast at 4.3 K the I-V curves 
can be fitted by the Nordheim- Fowler equation typical for 
tunnel emission, indicating a different transport mecha- 
nism at low temperatures Simple exponential fits 
and fits by power law were not appropriate to describe 
our data. From the measurements at different temper- 
atures we can deduce a negative temperature coefficient 
of resistivity which is indicative of thermally activated 
transport. Examining the differential conductance ([3]) 
gives a deeper insight into the charge transport mecha- 
nism. For T < AOK the electrons are tunneling when the 
applied field strength is sufficient to overcome the par- 
ticles charging energy. For T > AOK the electrons are 
thermally excited and thermally activated hopping can 
be expected. In the literature different hopping mech- 
anisms have been proposed: the simple thermally acti- 
vated hopping Q and variable range hopping (VRH) 
suggested for granular metals which are typical highly 
disordered systems 43]. For thermally activated trans- 
port the conductivity a scales with 

Here i?a is the activation energy and fee is the Boltzmann 
constant. For the VRH driven transport a different cor- 
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FIG. 5: Calculated capacitance values of a sphere in a hexago- 
nal lattice for different spheres radii and intersphere distances 
with corresponding numerical fits. The inset of the upper fig- 
ure shows the geometry used for FEM calculations, the hor- 
izontal sections of the calculated potential distribution, and 
the normalized electric field. 

relation was proposed [i^ : 

CTocexp^-^^ ■ (3) 

Here Tg can be interpreted as the activation tempera- 
ture. While the first mechanism describes the charge 
transfer between spatially adjacent particles, the VRH 
mechanism considers the transfer between the energeti- 
cally closest particles states. In order to distinguish be- 
tween these two cases we plotted the zero conductance 
data in an appropriate way for temperatures where the 
Coulomb blockade was overcome (see Supporting Infor- 
mation) . From the linear fits we conclude that the simple 
thermally activated transport model describes our data 
better than any other applied models. From the slope of 
the Arrhenius plot a value of 20 meV for the activation 
energy was obtained for the 7.6 nm particles. For smaller 
particles a value of 36meV was obtained while for the 
10.4 nm particles a value of 12meV was determined. In 
the literature [igj this activation energy was identified as 



the particles charging energy introduced in Equation [1] 
however, the Coulomb charging energy from Equation [1] 
should depend linearly on the particles' array width along 
the charge transport direction. Our data suggests instead 
that the activation energy scales with the particles' size. 
Small particles possess lower self-capacitance, which re- 
sults in higher activation energy according to Equation[T] 
In contrast, varying the interparticle distance did not af- 
fect the activation energy; our measurements always gave 
values in the range of 18-20 meV for the 7.6 nm par- 
ticles for the range of particle separations studied. To 
vary the interparticle distance we replaced the hexadecy- 
lamine ligands present after synthesis by dodecylamine, 
tetradecylamine and octadecylamine. For these cases, 
the interparticle distances as determined by GISAXS var- 
ied between 0.7nm and 1.6 nm. At room temperature the 
film conductivity was found to scale exponentially with 
the interparticle spacing (see 2]). 

The striking feature of the obtained data is the non- 
linearity of the TV curves even at higher voltages (see 
[3]). This observation encouraged us to verify the voltage 
dependence of the activation energy. In contrast to the 
results published so far we observed exponential decay of 
the activation energy at rising field strength (|4|). In fact, 
we find that the activation energy is not equal the charg- 
ing energy. At low fields electrons have to overcome the 
tunneling barrier and the particles' charging energy. Ris- 
ing field strength increases the electron energy until the 
Fermi energy of the electron sea in particle matches the 
energy of empty states of a nearby particle, into which the 
electrons are injected. In this regime the charging energy 
Ec is overcome by tunneling electrons, thermally assisted 
tunneling is not affected by further increase in voltage. 
From this point of view the particles charing energy can 
be calculated as the difference between the experimental 
zero field activation energy Eg and the high field limit 
activation energy Ei^f. Then the measured activation en- 
ergy can be expressed by the following equation. 

E, - Em + {Eo - Em) exp (^--^ ^ (4) 

^EM+{E,)exp(^-^^ (5) 

Here Uq represents a fitting parameter with voltage units. 
For small particles the Coulomb blockade was found to 
be 25.4 meV, for midsized particles 14.1meV, and for big 
particles 6.7 meV. 

Equation [T] suggests that one can calculate the charg- 
ing energy from a knowledge of Cs, the total capaci- 
tance "felt" by the tunneling electrons. Cs can be calcu- 
lated from the sum of the particles' self-capacitance Csoif 
and the interparticle junction capacitance Cgeo. The self- 
capacitance of a free standing sphere can be calculated 
analytically using equation [6] 

Cscif = ATTeoErr (6) 

Here Sr is the relative permittivity, Sq the vacuum per- 
mittivity of the interparticle environment and r is the 
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sphere radius. The junction capacitance depends on the 
geometry of the system under investigation and usuaUy 
cannot be calculated analytically. In the literature the 
model of a metallic sphere with dielectric shell embed- 
ded in metallic continuum has been used in many cases 
[H [Sli 113, SI after the paper of Abeles et al. 0]. Al- 
ternatively, the capacitance of a hemisphere over a plane 
surface has been approximately calculated [l^. These 
models deviate considerably from the experimental ge- 
ometry and potential distribution. In order to obtain the 
value of the Coulomb blockade we computed the geomet- 
rical capacitance by means of the finite element method 
(FEM). As the model a hexagonal lattice of spheres with 
7.6 nm diameter was used ([5]). This diameter corresponds 
to the medium-sized nanoparticles used in our experi- 
ments. A linear potential gradient was applied to the 
spheres. For the underlyin g ox ide layer a dielectric con- 
stant £r — 4.48 was used [431 which had only a slight 
influence on the obtained values. For the spheres envi- 
ronment the vacuum permittivity was used. The corre- 
sponding capacitance values for diverse dielectric envi- 
ronments can be calculated easily as the capacitance is 
proportional to the dielectric constant. Also a model with 
30° rotated potential gradient was used giving negligible 
deviations. The calculations gave a very low capacitance 
of 1.1 ■ 1Q~'^^]F. This result shows that in contrast to the 
results reported previously using less appropriate models 
the tunneling junction capacitance can be neglected com- 
pared to the self-capacitance of particles as was assumed 
previously by Averin and Korotkov in [soj ] . Nevertheless 
it is not possible to obtain the activation energy from 
the calculated capacitance as the value of the dielectric 
constant of the interparticle environment in the experi- 
mental setup is not known. 

In order to verify at which point the geometrical ca- 
pacitance starts to dominate over the self-capacitance, 
we calculated geometrical capacitances for different parti- 
cle sizes and physically reasonable interparticle distances. 
The resulting values of capacitance versus radius can be 
fit well by a parabolic equation, while the values obtained 
from the interparticle distance variation were fitted by a 
polynomial equation: 



C{S)^a + bS~^+c5-^. (7) 



In equation [7] a, 6, and c are fit parameters, and 5 
is the shortest distance between the particles' surfaces. 
From the calculated fit parameters we were able to ob- 
tain a universal equation which allows one to calculate 
the capacitance of a sphere in an hexagonal lattice with 
various sphere sizes and intersphere distances. If the rel- 
ative dielectric constant of the interparticle environment 
is known, the calculated capacitance has to be multiplied 



additionally by this value. 

Cgeoir, 6) = 6.25 • lO^^m-^Er {a + bS'^ + cS''^) (8) 

a = 8.6 • IQ-^^ F (9) 

6= -2.5- 10"^2 Fm (10) 

c 1.5 • 10"^^° Fm^ (11) 

For interparticle distances smaller than 1 nm, the spheres 
diameter has to exceed approximately IpLvn in order for 
the geometric capacitance to overcome the particles self- 
capacitance. 

In the present work we demonstrated the utilization 
of the Langmuir-Blodgett technique for deposition of 
cobalt-platinum nanoparticles of different sizes (4.4 nm, 
7.6 nm and 10.4 nm diameter) on structured silicon ox- 
ide substrates. The obtained ordered particle monolayers 
were characterized by SEM and GISAXS. 

The DC measurements on the films at temperatures 
between 50 and 300 K showed a simple thermally acti- 
vated electron hopping between spatially adjacent parti- 
cles with activation energies in the range of 12-36 meV. 
Small particles showed the highest values and big parti- 
cles the smallest. Calculations and experimental results 
suggest the charging energy to depend only on the par- 
ticles' self-capacitance and to be independent of the in- 
terparticle distance and thus independent of the geomet- 
rical capacitance, which was shown by FEM calculations 
to be negligible in the nanoparticle arrays. The geomet- 
rical capacitance starts to be predominant in comparison 
to the self-capacitance when the particles' diameter ex- 
ceeds one micrometer. Therefore the charging energy of 
smaller particles depends only on the dielectric constant 
of the interparticle medium and on the particles' radius. 

The I-V curves were fit by equations derived for semi- 
conductors, indicating the similarity between the trans- 
port mechanisms in semiconductors and metal nanopar- 
ticle arrays. Indeed the charging energy in the metal 
nanoparticle arrays plays an analogous role as the band 
gap in semiconductors. 

The investigated nanoparticle films showed an expo- 
nential decrease of the activation energy with the applied 
electric field. This dependence allowed us to calculate 
the particles charging energy, giving values in the range 
of 6.7 — 25AmeV. To our knowledge this approach has 
not been previously reported. 

The long range ordering of particles in our films gave 
rise to the observed effects which were explained by a 
simple next neighbor electron hopping. The size depen- 
dent effects are supported by numerical calculations. 
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